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Recent numerical and analytic computations based on the self-force (SF) formalism in general 
relativity showed that half-integral post-Newtonian (PN) terms, i.e. terms involving odd 
powers of 1/c, arise in the redshift factor of small mass-ratio black-hole binaries on exact 
circular orbits. Although those contributions might seem puzzling at first sight for conservative 
systems that are invariant under time-reversal, they are in fact associated with the so-called 
non-linear tail-of-tail effect. We shall describe here how the next-to-next-to-leading order 
contributions beyond the first half-integral 5.5PN conservative effect (i.e. up to order 7.5PN 
included) have been obtained by means of the standard PN formalism applied to binary 
systems of point-like objects. The resulting redshift factor in the small mass-ratio limit fully 
agrees with that of the SF approach. 


1 Introduction 

Stellar-mass compact objects inspiraling gradually about super-massive black holes may produce 
gravitational waves detectable by future space missions such as eLISA^. These systems, referred 
to as Extreme Mass Ratio Inspirals (EMRIs), can probe the strong gravity field regime, but 
proper data analysis of the resulting signal will require accurate waveform templates built from 
theoretical models. This has motivated, over the past ten years, numerous studies on the 
dynamics of point-like objects on a curved background Due to the metric perturbations 

generated by its own mass-energy, the point mass effectively feels a self-force (SF) that induces 
deviations from the geodesic wordline followed by a test particle on the background. When 
the background is a black-hole spacetime, the acceleration may be sought in the form of an 
asymptotic expansion in powers of the mass ratio q = mil m 2 1. 

In the SF approach, the first order perturbation of the background metric is 
obtained by convolving over the stress-energy tensor the regularized (R) Green function 
{x,x') that solves the linearized homogeneous Einstein equations in harmonic gauge and 
has the property that it coincides with the retarded Green function when x lies in the chrono¬ 
logical future of x', while vanishing when x is in the chronological past of x' The trajectory 
of the particle is then precisely that of a geodesic for the perturbed metric + 5g^y. 

By contrast, the post-Newtonian (PN) approach is based on the formal expansion, on a flat 
background spacetime, of all quantities of interest, in powers of u/c, where v is the largest typical 
velocity of the problem. The standard PN approach is first defined for general extended PN 
sources with compact support and then specialized to compact binary systems®. In that case. 


V is taken to be the relative coordinate velocity vi2- Moreover, when the bodies are compact, 
they may be effectively represented as point particles. Ultra-violet (UV) divergences at the 
particle positions are tackled by means of dimensional regularization The PN expressions 
are valid at a given coordinate time in a spatial region, referred to as the near zone, that entirely 
contains the matter source and whose radius is much smaller than the gravitational wavelength. 
Because SF and PN methods are so radically different from each other, notably regarding the 
regularization schemes, comparing PN expansions of observable quantities truncated at linear 
order in q to their SF counterparts expanded in power of 1/c allows for non-trivial cross-checks 
that strengthen our confidence in both perturbative techniques. 

After the first comparison between PN and gravitational SF calculations^, rapid progress has 
been made over the last six years, mainly due to both high precision numerical computations from 
a SF perspective and extensive analytical PN computations Recently, after the possibility 

for this comparison had been dramatically extended from the SF side it was realized that 
observable quantities could contain half-integral ^ PN terms that are nevertheless conservative, 
starting at the order 5.5PN . Here, we shall explain how those terms can arise within the PN 

framework and sketch their calculation at the next-to-next-to-leading order . As we shall 
see, they are closely related to the so-called tail-of-tail effects in general relativity. The success 
of this SF/PN comparison actually provides an excellent test of the intricate PN machinery for 
computing non-linear tail-of-tail effects — these being relevant for template waveform generation 
of comparable mass compact binaries to be analyzed in ground and space based detectors. 


2 Comparing post-Newtonian and self-force results 


2.1 The Detweiler variable 

We shall focus on the Detweiler variable^, which represents physically the inverse of the redshift 
of a photon emitted by a particle moving on an exact circular orbit around a Schwarzschild 
black hole and detected by an infinitely far-away observer along the rotation axis. The ensuing 
spacetime is helically symmetric, with a helical Killing vector K" tangent to the four-velocity 
uf on the particle worldline. Alternatively, the redshift variable is defined geometrically as the 
conserved quantity associated with the helical Killing symmetry relevant to spacetimes with 
exactly circular orbits. In an appropriate class of coordinate systems, the redshift factor uj 
reduces to the t component of the particle’s four-velocity. 
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where gafsiui) denotes the metric evaluated at the particle’s location yf = {ct,y\) by means of 
dimensional regularization, and vf = dyf/dt = {c,v\) is the coordinate velocity. 

The Detweiler variable (1) has been computed to high-order using on the one hand standard 
PN theory suplemented with dimensional regularization, valid in weak field and on the 

other hand both numerical and analytical SF approaches, valid in the limit q = mijm 2 1. 
Over the last two years, its accuracy has improved drastically on the SF side due to the new 
application of methods to represent analytic solutions for metric perturbations of black-hole 
spacetimes. In a first stage, based on some exact solutions of the Teukolsky equation the 
PN coefficients of the redshift factor were obtained numerically to 10.5PN order Analytic 
expressions were even found for a subset of coefficients, specifically those that are either rational, 
or made of the product of vr with a rational, or a sum of commonly occurring transcendentals. 
An alternative SF approach^®, based on the post-Minkowskian expansion of the Regge-Wheeler- 
Zerilli equation ', has also reached PN coefficients analytically up to 8.5PN order. Most recently, 
both methods have been extended to extremely high orders, typically 21.5PN for the redshift 
factor 




The appearance of half-integral PN coefficients (of type ^PN with n being an odd integer) 
in the conservative dynamics of two particles on circular orbits is a feature of high-order PN 
expansion. Resorting to standard PN methods, we shall show now that the leading half-integral 
PN terms originate from non-linear integrals depending on the past history of the source — 
so-called hereditary type integrals. 

2.2 Dimensionality argument 

The fact that terms at half-integral PN orders cannot stem from the source variables evaluated 
at the current time follows from the general structure of “instantaneous” terms entering the 
redshift factor (1) in the center-of-mass frame. After replacing all accelerations by the lower- 
order equations of motion, takes the form (with usual Euclidean notation) 



integers 


where m = mi m 2 , n = mim 2 /m^, while ni 2 = {yi — y 2 )lri 2 stands for the relative direction 
of the two particles. We have taken the expansion when the mass ratio —)• 0. 

The counting of the 1/c powers shows that the PN order of the generic term in Eq. (2) can be 
half-integral only if p is odd, in which case it vanishes for circular orbits, when the velocity V 12 
and unit direction ni 2 are evaluated at the same current time t. However, integration over some 
intermediate time extending from the infinite past up to t could allow a coupling between these 
vectors at different times. Of course, in general relativity, this type of “hereditary” dependence 
over the past of the system does occur, due, in particular, to wave tails produced by the back- 
scatter of linear waves on the spacetime curvature. 

3 Post-Newtonian compntation of half-integral order contributions 

3.1 Structure of the tail contributions 

The tail effects, associated with non-linear wave propagation, are best investigated by construct¬ 
ing the (multipolar-)post-Minkowskian expansion of the metric in powers of the gravitational 
constant G, outside the matter source®. We start from the Einstein field equations in vacuum, 
written in terms of the field perturbation variable h^'^ = y/—gg^'^ — on the flat background 
77 ^*^ in Cartesian coordinates {x*}, with g representing the determinant of g^p. Adopting the 
harmonic-gauge condition dyh^'^ = 0, the relaxed Einstein equations for reduce to the wave¬ 
like equations = A^^, with □ = y'^^dadp. The non-linear source term is an expression 

of second-order (at least) in At linear order in G, the most general solution depends on six 
sets of source multipole moments: the mass-type moments, {i being the multipole or¬ 

der), the current-type moments Jl = Tq,,.i£, and four sets of so-called gauge moments, irrelevant 
for the present discussion. The higher order solutions are obtained iteratively by applying the 
flat retarded integral operator on the source term, after multiplication by a regularization 
factor r^ to cope with the divergence of the multipole expansion when r = |x| —)• 0 . Analytic 
continuation in R G C is invoked and the finite part (FP) when R —)• 0 provides a certain 
particular solution. To ensure that the harmonic coordinate condition is satished at each step, 
we must add to the latter solution a specific homogeneous retarded solution, which does not 
generate tail integrals and can be safely ignored here. 

Since, ultimately, we shall be interested in the metric at the location of one of the particles, 
our goal is to compute the near-zone expansion, indicated below by an overline, of the general 
solution initially defined in the exterior of the source, when r —)• 0. It is obtained directly at a 
given order from the near-zone expansion of the corresponding source, without need to control 
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the full solution, from the formula 


J [hLS{r, t - r/c)] = di 
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where hi denotes the symmetric trace-free part of ul = {£ G N). The first term is a 

homogeneous solution of the wave equation which is of retarded-minus-advanced type and regular 
at r = 0. It may be directly expanded in the near-zone, where it is valid by virtue of a matching 
argument. The second term in Eq. (3) is a particular solution of the inhomogeneous equation 
that is defined by means of the “instantaneous” inverse box operator, representing the 

PN expansion of the symmetric integral operator supplemented with a regulator r^ multiplying 
the source and a finite part as i? —)• 0. This term diverges when r —)• 0 and should be matched 
to a full-hedge solution of the held equations inside the source. However, we proved (see the 
appendix there) that it cannot actually contribute at half-integral PN orders, so that the ehect 
we are looking for comes only from the term containing ^(u). The latter function is given by a 
specihc double integral over the source piece S(r, t), and regularized by the hnite part as H —)• 0, 


where R{p, s) = f dA S{X, s). (5) 

B Jo 

This function is the crucial object to investigate for the purpose here. 

At linear order in the mass ratio, we may disregard any hereditary term involving the 
product of more than two moments other than the mass monopole M, since each such multipole 
is proportional to v. General results on the structure of the gravitational field in harmonic gauge 
tell us that hereditary contributions of type M x • • • x M x A4p, with J\Ap = Ip or Jp, read® 
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where the upper sign (a) refers to time derivatives and the tensor function K^p{t, u) is a dimen¬ 
sionless kernel. Using dimensional analysis combined with “angular-momentum” selection rules, 
it is straightforward to show that only interactions with k ^ 3 can produce the half-integral PN 
terms of interest, starting at the leading 5.5PN order. In fact, at the next-to-next-to-leading 
order beyond 5.5PN, we can restrict ourselves to hereditary cubic interactions M x M x JAp, 
which may be interpreted physically as gravitational-wave tails of tails. They must be computed 
at the 5.5PN, 6.5PN and 7.5PN orders for the mass qnadrupole, 6.5PN and 7.5PN orders for 
the mass octupole and current quadrupole, and so on. The leading contributions of current-type 
moments are IPN order higher than those of the mass moments. 


3.2 Sketch of the calculation of uj 

The elementary source terms hLS{r,t — r/c) for the tails of tails can be either instantaneous, 
with S{r,t—r/c) = r^JAp\t—r/c), or hereditary. In the latter case, S{r,t—r/c) is an integral of 
the type r~^ Qmix)JAp'^ [t — rx/c), with Qm{x) being a Legendre function of the second 

kind. After some transformations, the tail-of-tail piece of Q{u) can be written as the finite part at 
H = 0 of some coefficient C^km{B), times an integral of the variable r whose integrand involves 
the regulator times derivatives of JAp{t — r). Now, we find that Ci^rniB) may comprise 
(simple) poles at the order we are working, so that the factor generates a logarithm kernel 
Inr. Insertion of the former piece of Q into the homogeneous wave entering Eq. (3) yields the 
form (6), with K^p{t,u) proportional to ln(t — u) for t > u, and zero elsewhere 









At this stage, it is important to realize that this sort of “pure” tail-of-tail contributions can 
generate another kind of half-integral terms, “mixed” contributions, by coupling to the non¬ 
tail part of the PN metric in the source of Einstein’s equations. This part is obtained most 
conveniently by solving the relaxed field equations in the near zone, where is augmented by 
the matter source 167rGc“'^|g(|r^'^. It is usually parametrized by means of appropriate potentials, 
such as the Newtonian potential U = Gmijri + Gm 2 jr 2 (with ri = |a; — yi|). The time 
component of the gravitational field, for instance, is composed of “ordinary” PN terms: 
—AU/(P‘ -!-•••, plus tail terms containing our effect: -!-•••. Its product with, say, 

whose structure is similar, produces couplings that must be crucially taken into account 
in the calculation. Their number is minimized by moving to an adapted gauge. Quadratic 
and cubic PN iterations are then required to find the complete half-integral PN part of the 
metric at the 7.5PN order. The successive solutions are constructed by means of hierarchies 
of “superpotentials” derived from the potentials that enter at the 2PN order In the 
end, we decompose the tail integrals into conservative time symmetric and dissipative time 
anti-symmetric pieces and simply discard the dissipative piece from our results. 

Using the standard stress-energy tensor for point particles to model the binary, all multipole 
moments, potentials and superpotentials can be evaluated explicitly for circular orbits. The 
hereditary integrals are derived from standard formulas. This yields for the redshift 
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where we have written only the relative 2PN terms relevant to our next-to-next-to-leading order 
calculation, with all the other terms indicated by ellipsis. The result (7) is in full agreement with 
those derived by gravitational SF methods, either semi-analytical or purely analyticaU^- We 
emphasize that it has been achieved from the standard PN approach, which is not tuned to a 
particular type of source (contrary to various analytical and numerical SF calculations), as it is 
actually applicable to any extended PN source with compact support. The 7.5PN order reached 
by the present calculation is perhaps the highest order ever reached by traditional PN methods. 
Note also that while SF results may be relativistic but with q <C 1, the present method, valid 
primarily in the PN regime, is in principle applicable for arbitrary mass ratios. The time is now 
ripe for the SF approach to proceed to second order in q. We look forward to an occasion when 
high precision SF methods applied at second order may be fruitfully compared with known PN 
results in the weak field regime. 
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